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ABSTRACT 

We present a 9 deg^ map of the North American and Pehcan Nebulae regions obtained in all four 
IRAC channels with the Spitzer Space Telescope. The resulting photometry is merged with that at 
JHKg from 2MASS and a more spatially limited BVI survey from previous ground-based work. We 
use a mixture of color- color diagrams to select a minimally contaminated set of more than 1600 objects 
that we claim are young stellar objects (YSOs) associated with the star forming region. Because our 
selection technique uses IR excess as a requirement, our sample is strongly biased against inclusion of 
Class HI YSOs. The distribution of IRAC spectral slopes for our YSOs indicates that most of these 
objects are Class II, with a peak towards steeper spectral slopes but a substantial contribution from 
a tail of flat spectrum and Class I type objects. By studying the small fraction of the sample that is 
optically visible, we infer a typical age of a few Myr for the low mass population. The young stars 
are clustered, with about a third of them located in eight clusters that are located within or near the 
LDN 935 dark cloud. Half of the YSOs are located in regions with surface densities higher than 1000 
YSOs / deg^. The Class I objects are more clustered than the Class II stars. 

Subject headings: stars: formation - stars: circumstellar matter - stars: prc-main sequence - ISM: 
clouds - ISM: individual (NGC7000, IC5070) - infrared: stars ~ infrared: ISM 



1. INTRODUCTION 

Much of our current knowledge regarding star-forming 
patterns and circumstellar disk evolution derives from 
study of molecular cloud complexes within a few hun- 
dred parsecs of the Sun. Among these are a large num- 
ber of lower-mass clouds such as Taurus and, more in- 
frequently, dense clouds like Orion, which is the proto- 
typical high-mass and high density star forming region. 
While nearby cloud complexes serve as our primary em- 
pirical guide to understanding the formation and early 
evolution of stars, it is important that we study more 
than just the nearest examples. One little-studied-to- 
date star-forming region beyond the solar neighborhood 
is the North American and Pelican Nebulae (NGC 7000 
and IC5070, respectively), towards Z=85, 6=— 0.5, about 
600 pc from Earth. This region is probably more repre- 
sentative of star formation in the disk of the Milky Way 
than Taurus or Orion, it is also the next closest region 
after Orion with a substantial numbers of intermediate 
and high mass stars. The physical appearance of the 
nebulae on optical images is thought to reflect a combi- 
nation of a large, background H II region (W80) overlaid 
by several foreground, dark clouds, with the edges of the 
dark clouds illuminated by the optically-unseen primary 
exciting star of the H II region. Because these nebulae lie 
essentially in the plane of the Galaxy and, as seen by us, 
along a spiral arm, study of the region is further confused 
by the juxtaposition of young stars and star-forming re- 
gions at a variety of distances. For the purposes of the 
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rest of the paper, wc will refer to the entire region as the 
N ANeb, o r the NAN complex. 

iHerbid (fl95l was one of the first to study the star 
formation associated with these nebulae. He identified 
68 Ha emission line stars from objective-prism plates 
of the region, including a small cluster of pre-main se- 
quence (PMS) stars in the "Gulf of Mexico" portion of 
the NANeb. Many of the young stars identified by Her- 
big lie along the bright rims of the dark clouds, possibly 
indicating that they are recent, second generation star 
formation products triggered by the expanding H II re- 
gion from the first generation of stars in the NANeb. 
Herbig estimated a rough distance of 500 pc for the star- 
forming complex, but with a variety of caveats related to 
whether this distance applied to some or all of the stars 
and nebulae. The best avail able modern distance for th e 
region appears to be that of iLaugalvs fc Straizvsl (|2002f l. 
who estimate a distance of ~600 pc. We will use this 
distance for the purposes of the paper with the caveats 
n oted above. 

iBallv fc Scovilld (|1980f ) mapped the CO associated 
with the dark clouds that define the Gulf of Mexico and 
Atlantic Ocean, and estimated that the mass in gas still 
present is of order 3 x 10^ Mq. Their interpretation of 
the CO data was that the remaining dark clouds are the 
remnants of a much larger molecular cloud complex that 
has largely been disrupted by the hot, massive stars that 
were f ormed a few million years ago. iCambresv et al.l 
(I2OOI used the Two -Micron All-Sky Survey (2MASS; 
iSkrutskie et al.l 120061 ) point-source catalog to map the 
extinction and star clustering towards this region, deriv- 
ing extinctions of > 35 magnitudes in some portions 
of the darkest cloud (Lynds935). Cambresy et al. also 
i dentified nine appar e nt sta r clusters in the 2MASS data. 

IBallv fc ReipurthI (|2003f ) obtained narrow-band Ha 
and [S II] imaging of the region of the Pelican nebula, and 
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TABLE 1 

Summary of IRAC observations (programs 20015 and 462) 



Program : field-ID 



map center 



AORKEY 



P20015 


12 


P20015 


13 


P20015 


21 


P20015 


22 


P20015 


23 


P20015 


24 


P20015 


31 


P20015 


32 


P20015 


33 


P20015 


34 


P20015 


42 


P20015 


43 


P00462 





P00462 


1 


P00462 


2 


P00462 


3 


P00462 


4 


P00462 


5 


P00462 


6 


P00462 


7 



20h55m30.00s, 
20h53m06.00s. 
20h57m54.00s, 
20h55m30.00s, 
20h53m06.00s. 
20h50m42.00s. 
20h57m54.00s, 
20h55m30.00s. 
20h53m06.00s, 
20h50m42.00s. 
20h55m30.00s. 
20h53m06.00s, 
20h50m38.00s. 
21h01ml3.00s. 
20h58m51.00s, 
20h59m05.00s. 
21h00ml7.00s. 
21h01m50.00s, 
20h50ml9.00s. 
20h48m44.00s, 



+44d59m00.0s 
+44d59m00.0s 
+44d23m00.0s 
+44d23m00.0s 
+44d23m00.0s 
+44d23m00.0s 
+43d47m00.0s 
+43d47m00.0s 
+43d47m00.0s 
+43d47m00.0s 
+43dllm00.0s 
+43dllm00.0s 
+45d05ml0.0s 
+44d54m43.0s 
+45dl2m42.0s 
+42d59m38.0s 
+43d23m25.0s 
+44dl0m34.0s 
+42d58m37.0s 
+43dl7m51.0s 



16790528 
16790784 
16791040 
16791296 
16791552 
16791808 
16792064 
16792320 
16792576 
16792832 
16793088 
16793344 
24251392 
24251648 
24251904 
24252160 
24252416 
24252672 
24252928 
24253184 



identified a number of new HH objects and coUimated 
outflows, mostly originating from the interface region be- 
tween tlie dark cloud and the H II region. In a couple of 
cases, the flows originate from near the tip of elephant 
trunks. The exciting sources for most of the flows could 
not be identified from the existing data. lArmond et al.l 
(|2003f ) obtained additional narrow and broad band op- 
tical and near-IR imaging of the NANeb, and identified 
28 more HH objects, including many associated with the 
PMS stars of Herbig's Gulf of Mexico cluster. More re- 
cently, the probable exciting source for the H II region 
has been identified as an 05 star behind about Av=lO 
mag of extinction, located inside the dark cloud LDN 935 
separating the North Amer ica and the Pelican Nebulae 
()Comer6n fc Pasauali|[2005l ). 

We have conducted a large (~9 deg^) i nfrared imag- 
ing s urvey with Spitzer Space Telescope (| Werner et al.l 
[200l of this region, along with supporting data obtained 
in the optical for the ^ 2.4° x 1.7° central region. This 
present paper, the first of a series, presents the Infrared 
Array Camera (IRAC) data. Future papers will cover 
the Multiband Imaging Photometer for Spitzer (MIPS) 
data, with special emphasis on the Gulf of Mexico clus- 
ter (RcbuU et al. 2009; hereafter R09), and the optical 
classification spectroscopy (Hillenbrand et al. in prep). 

First, we present the observational details for the IRAC 
observations (ij2]). We then use the IRAC colors to se- 
lect a minimally contaminated, if not complete, sam- 
ple of YSO candidates and discuss their proper- 
ties (classes, spatial distribution) in the context of other 
star-forming regions. Because IRAC is very sensitive to 
emission from HH objects, we then discuss the IRAC ob- 
servations of one of the HH objects found in this region 

m- 

2. OBSERVATIONS, ANCILLARY DATA, AND BASIC DATA 
REDUCTION 

2.1. Observations 

The Spitzer observations of the NAN complex were 
obtained as part of the Cycle-2 program 20015 (PI: L. 
RebuU). Additional observations of the "corners" of the 



map were obtained as part of program 462 in an effort 
to increase the legacy value of the data set. The IRAC 
observations from program 20015 were obtained 9-11 Au- 
gust 2006, amd the IRAC observations from program 462 
were obtained 15-27 November 2007. IRAC (jFazio et al.l 
120041 ) observes at 3.6, 4.5, 5.8, and 8 ^m. 

The Cycle-2 IRAC observations were designed to cover 
the region of highest extinction in a manner as indepen- 
dent of observing constraints as possible. Figure [T] is the 
NANeb region in the optical (from the Digitized Palo- 
mar Optical Sky Survey). The region covered by our 
IR AC map is indica t ed, as is the = 5 contour from 
the iCambresv et al.l (|2002f ) extinction map. The obser- 
vations were broken into 12 astronomical observation re- 
quests (AORs); the AORKEYs are given in Table[TJ This 
region is at a high ecliptic latitude (+57°), so the field of 
view rotates by about a degree a day, necessitating small 
individual AOR coverage in order to completely tile the 
region without leaving gaps. Also, because of the high 
ecliptic latitude, asteroids are not a significant concern, 
and we therefore obtained all of the imaging at a single 
epoch. Each mapping AOR was constructed with the 
same strategy - at each map step, four dither positions 
were observed, each with high-dynamic-rangc (HDR) ex- 
posures of 0.6 and 12.0 second frame times. For this 
strategy, the on-line SENS-PET for a high-background 
region reports 3-cr point-source sensitivities of 7.2, 10.7, 
66, and 78/xJy for the 4 IRAC bands, respectively. Dif- 
ferential source count histograms for our point sources 
are very similar to those for other star-forming regions 
observed in a similar manner - see in particular Figure 8 
of Winston et al. 2007 (ApJ 669, 493) for 3.6 iim; the 
NaNeb histograms indicate that our 90% completeness 
limits are about [3.6] = 15 decreasing to about [8.0] = 
12.5. For our purposes, however, the more important 
number is the completeness for detecting objects in all 
four channels. We estimate this by determining where 
the four-channel differential source counts drops below 
the 3.6 source counts by more than 20%, which oc- 
curs at about [3.6] = 12.2. This is about as expected 
given that four-channel catalog is normally limited by 
the detections in the 8.0 fim channel. This completeness 
limit corresponds to a spectral t ype of M4 or Ml acc ord- 
ing to the BCAH98 isochrones (jBaraffe et al.l 119981) for 
1 or 5 Myr, respectively, a distance of 600 pc (assuming 
no reddening) and the temperature to spectral type scale 
defined by Luhman et al. (2003). The faintest object in 
our final catalog has [3.6]^16, corresponding to a mass 
of 0.02-0.03 M0 at 1-5 Myr age, assuming the 3.6 fim flux 
is photospheric. 

The remaining IRAC observations from program 462 
were designed using the fixed-cluster observing mode to 
cover the irregularly-shaped regions to the same depth 
as the main map and to create a final map that is ap- 
proximately square in shape. The map center given in 
Table [T] is the approximate center of each AOR. We used 
the software developed by R. Gutermuth and T. Megeath 
(private communication) for use by the IRAC GTO and 
Gould's Belt teams to construct these AORs. 

2.2. Basic Data Reduction 

We started with the Spitzer Science Center (SSC) 
pipeline-produced basic calibrated data (BCDs), version 
S14.4. We ran the IRAC Artifact Mitigation code written 
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Fig. 1. — Optical image (~ 4 X 3.2 degrees) of the NANcb region from the Palomar all-sky survey. North is up, East is to the left. The 
left and right edges arc at 21h02m32s and 20h46m52s; the top and bottom edges are a t 45d25m40s and 42d40 m56s. The region observed by 
IRAC is indicated by red lines. The green line is a smoothed Av=5 contour from the lCambresv et al.l I I2002I ) extinction map. The galactic 
coordinate of the center of this image is 1=84.8, b=-0.6. 



by S. Carey and available on the SSC website. We con- 
structed mosaics from the corrected BCDs using the SSC 
mosaickin g and point-source extracti on software package 
MOPEX (jMakovoz fc Marleaul[2005l ). The mosaics have 
a pixel scale of 1.22" px~^, very close to the native pixel 
scale. Figure [2] shows the mosaics in channel 1 (3.6 ^m), 
and Figure [3] shows a 3-color IRAC mosaic (4.5, 5.8, and 
8 /im). 

We performed aperture photometry on the combined 
long- and short-exposure mosaics separately using the 
output of the APEX (par of the MOPEX software) de- 
tection algorithm (APEX is a part of the MOPEX pack- 
age) and the "aper" IDL procedure from the IDLASTRO 
library. We used a 2 pixel radius aperture and a sky an- 
nulus of 2-6 pixels. The (multiplicative) aperture correc- 
tions we used follow the values given in the IRAC Data 
Handbook: 1.213, 1.234, 1.379, 1.584 for IRAC channels 
1, 2, 3, and 4, respectively. Fluxes have been converted to 
magnitude given the zero magnitude fluxes of 280.9 ±4.1, 



179.7±2.6, 115. Oil. 7 and 64.1 ±0.9 for channels 1, 2, 3, 
and 4 respectively (IRAC Data Handbook). We take our 
photometry errors to be that given by the IDL proce- 
dure. We have checked that the uncertainties derived by 
the IDL procedure are consistent with the dispersion in 
IRAC colors for bright, off-cloud sources. For the longer 
integration times, the estimated average errors on our 
photometry in the 4 IRAC channels are 0.020, 0.025, 
0.031, 0.035 magnitudes for bright stars, increasing to 
--0.1 mag for channels 1-4 at 15, 14.7, 13.2, 12.7 mag. 

Since the sources have been extracted individually 
from the mosaiced BCDs, we have used the overlap be- 
tween BCDs to estimate our astrometric precision. The 
histogram of coordinate differences shows a one sigma 
RMS uncertainty of 0.3" (~ 1/4 of pixel) in both direc- 
tions. 

The APEX source detection algorithm has a tendency 
to identify multiple sources within the PSF of a single 
bright source which can cause significant confusion at 
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the bandmerging stage. Because our nominal fluxes are 
from the 2 pixel aperture photometry, any object which 
has a companion within 2 pixels will have a confused flux. 
The photometry lists were therefore cleaned of multiple 
sources prior to the bandmerging. Objects with a com- 
panion within 2 pixels had one of the sources removed. 
The source chosen for removal was the one with the low- 
est signal-to-noise ratio. 

We extracted photometry from the long and short ex- 
posure mosaics separately for each channel, and merged 
these source lists together by position using a search ra- 
dius of 2 pixels (2.44") to obtain a catalog for each chan- 
nel. The magnitude cutoff where we transition to using 
the long rather than the short exposure photometry cor- 
responds to magnitudes of 11, 10, 8.4, and 7.5 magni- 
tudes for channels 1, 2, 3 and 4 respectively. 

Because of the complex nebulosity in this region, 
APEX (like any point-source detection algorithm) can be 
fooled by structure in the nebulosity. This effect is most 
apparent at 8 /xm, where there are also the fewest stel- 
lar (point) sources detected. To attempt to remove false 
sources, we have compared photometry obtained via a 2 
and 3-pixel aperture. After applying the aperture correc- 
tion given in the IRAC Data Handbook, we eliminated 
sources with a difference of magnitude > 0.2 mag; doing 
so rejected 32% of the raw detections in IRAC's 8 fim 
band. This process is only applied to the 8 /im channel 
because the nebulosity is most prominent there (due to 
the presence of strong PAH bands at 7.6 and 8.6 ^m). 
Also, because there are the fewest point sources detected 
in this band (see Table [5]), the chances of there being 
two legitimate sources within 2-3 pixels of each other are 
much lower than in channels 1 or 2. Spot-checking the 
images confirms this assertion. We acknowledge this step 
may eliminate some true YSOs from our 4-band catalog. 
We deem this to be acceptable because our goal is to 
produce a minimally contaminated catalog of YSOs, not 
a complete catalog. 

2.3. Ancillary Data & Bandmerging 

Table [2] provides statistics for our entire catalog of ob- 
jects; an object is included in our master catalog even if 
it is only detected in a single band. In order to create 
our final multi-wavelength catalog which we will use to 
identify new YSOs, we first merged the four individual 
IRAC source lists together, starting with IRAC-1, tak- 
ing the closest source within 1" as the best match. The 
radius of 1" has been chosen based on our astrometric 
precision of 0.3" (sec section above) and the density of 
sources in the 3.6 /xm image {i.e. the most crowded im- 
age). 

In order to provide photometry at other bands, we 
cross-matched our catalog to 2MASS, taking the closest 
source within 1", and then to our optical catalog, again 
taking the closest source within 1". Out of the 63 084 
objects detected at all 4 IRAC channels, 93% have a 
2MASS counterpart. For the entire region, 10% of the fi- 
nal catalog stars have an optical counterpart; out of just 
the region covered by the optical photometry data, 28% 
have an optical counterpart. The MIPS catalog will be 
discussed in detail in R09. In summary, we have covered 
approximately the same area as the IRAC map, and we 
performed the source extraction using APEX. 

BVI images were obtained with the KPNO 0.9 m tele- 



scope on four photometric nights in June 1997. An area 
of approximately 2.4xl.7deg^ was covered by mosaic- 
ing the 23.2' x 23.2' FOV CCD in a grid with overlap 
of typically 3' along each border. Exposures of 10 sec 
and 500 sec enabled unsaturated photometry between 
V = 11 and 20 mag. In IRAF, images were bias sub- 
tracted and flat fielded with sky fiats taken each night. 
Sources were identified using a 5a detection threshold 
and the DAOFIND task. The photometry was measured 
through an aperture of 5 pixels radius with background 
determined as the median in an annulus from 7-14 pixels 
where the plate scale was 0.68"/pixcl. Aperture correc- 
tions were applied. Absolute photometric calibration was 
achieved through observation of Landolt standards and 
further self-calibrated across the overlap regions using 
stars in common between frames. First the short frames 
were calibrated to the long frames then the long frames 
were tied across the distinct telescope pointings, ignoring 
individual photometric errors larger than 0.1 mag. The 
average of the 43 spatial frame-to- frame offsets was 0.024 
mag with dispersion 0.041 mag and a few frames requir- 
ing offsets as large as 0.1 mag. The individual offsets 
themselves have standard deviation typically 0.02-0.03 
mag and standard deviation in the mean <0.01 mag, 
which we take as the self-calibration error. Photome- 
try was adopted from the long frames except when in 
the nonlinear or saturated regimes of the CCD response 
in which case the short frame photometry was adopted; 
these criteria were applied at B <15.5, V <14.5, / <13.5. 
Astrometry was obtained using the HST Guide Star Cat- 
alog and the TASTROM task and is estimated accurate 
to <0.3" (Ict) based on stars in the overlap regions. 

Additional RI frames were obtained on a fifth night 
over a more limited area, 2x1.5 deg'^. These were reduced 
in a similar way, requiring a 5a detection threshold. 

3. IDENTIFICATION AND CHARACTERIZATION OF OUR 
YSO SAMPLE 

3.1. Selection of YSO candidates 

Several studies in the literature, usually focusing on 
subrcgions of this complex, have identified a total of 170 
YSOs in the area covered by our IRAC map, using a 
variety of techniques such as NIR excess or Ha emission. 
Most of these objects are earlier than K7. A complete 
list of the new Spitzer data for these previously-known 
YSOs will appear in R09. 

Now, with our new, comprehensive multi-wavelength 
view of the complex, we can begin to create a census of 
the members of this complex having infrared excesses. 
However, doing so is difficult, as it requires an exten- 
sive weeding-out of the galactic and extragalactic con- 
taminants. To begin this task, we have opted to cre- 
ate a minimally contaminated sample of Spitzer-selected 
YSO candidates, as opposed to identifying every possible 
YSO candidate. By miniminally contaminated sample, 
we simply mean a sample which includes as few non- 
YSOs as possible {i.e. from which AGB stars, AGN, 
and any other objects whose IR colors mimic YSOs have 
been eliminated) . We identify only stars with excesses at 
IRAC wavelengths as YSOs - hence any member without 
excess {i.e. Class HI YSOs) or stars with excess only at 
wavelengths greater than 8 fj,m. (inner cleared regions or 
transitional disks) will not be included. We now discuss 
the selection criteria we have used to select our reliable 
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Fig. 2. — Mosaic of the NANeb region from IRAC channel 1 (3.6 /^m). Orientation and map center are the same as previous Figure; 
overlays from previous Figure are included for reference. 



TABLE 2 

Numbers of IRAC sources extracted from the NANeb data. 



Band Number of sources 



Fractional number of sources wrhich match (in %)^ 
[3.6] [4.5] [5.8] [8.0] [3.6 + 4.5] [3.6 + 4.5 + 5.8] [3.6 + 4.5 + 5.8 + 8.0] 



3.6 

4.5 

5.8 

8.0 

JHKs 

BVI 

[24] 



558053 




75 


32 


12 


75 


30 


11 


510922 


82 




33 


13 


82 


33 


12 


200875 


90 


83 




32 


83 


83 


31 


77898 


84 


83 


82 




83 


81 


81 


220799 


90 


86 


65 


27 


86 


64 


26 


16595 


85 


83 


72 


40 


83 


71 


39 


4232 


51 


51 


48 


43 


51 


48 


43 



^ As an example of how to interpret this table, the "75" in cell ([3.6], [4.5]) means that 75% of the stars detected at 
3.6 f-LUi also are detected at 4.5 /-im 
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Fig. 3. — A 3-color view of the IRAC NANeb with 4.5 fim in blue, 5.8 fim in green, and 8 fim in red. Orientation and map center are 
the same as Figure [T] overlays from Figure [T] are included for reference. 
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YSO sample. 

Our selection method is described below and illus- 
trated in Figure ID We initially require detection at all 
four IRAC bands, which substantially limits the catalog 
(see Table [5]) and then we follow the IRAC four-band 
source characterization described in IGutermuth et all 
()2008l section 4.1). Using their extragalactic contam- 
ination criteria, we have first rejected 272 sources as 
having colors consistent with galaxies dominated by 
PAH emission. This selection has been made in the 
[4.5] - [5.8]/[5.8] - [8] and [3.6] - [5.8]/[4.5] - [5.8] planes, 
where we reject objects seen in the grey zones plotted in 
panel (a) and (b) of Figure [H this selection is defined by 
two sets of inclusive equations: 



[4.5] 



[5.8]< W([5-8]-[8.0]-l) 
& 

[4.5] - [5.8] < 1.05 

& 

[5.8] - [8.0] > 1 

OR (1) 



[3.6] - [5.8] < ^^^([4.5] - [8.0] - 1) 
& 

[3.6] - [5.8] < 1.5 
& 

[4.5] - [8.0] > 1 

We then rejected 823 more sources having colors con- 
sistent with AGN in the [4.5] vs. [4.5] -[8.0] plane; the 
region of color space used to select AGN-likc sources is 
plotted in panel (c) of Figure S] and is defined by the 
following equations: 

[4.5] - [8.0] > 0.5 
& 

[4.5] > 13.5 + ([4.5] - [8.0] - 2.3)/0.4 

& 

[4.5] > 13.5 



& 



[4.5] > 14+ ([4.5] - [8.0] -0.5) 
OR 

[4.5] > 14.5 - ([4.5] - [8.0] - 1.2)/0.3 
OR 
[4.5] > 14.5 



(2) 



(see Appendix A of lGutermuth et aI1l2008l for more de- 
tails). Note that the AGN-like source selection has 
been ma de in the observed [4.5] vs. [4.5] — [8.0] plane, 
wherea s IGutermuth et al. (2 0081) use the dereddcned 
plane. IGutermuth et alT ( 200'8h are able to work easily 



in the dereddened plane because they have a high spa- 
tial resolution map. Because the NAN complex is 
much further away than NGC 1333, we do not have the 
ability to construct such a high-resolution Ay map. We 
tested the implication s of this limitation by using the 
ICambresv et al.l (|2002( ) extinction map to deredden in- 
dividual objects. Just 25 sources would be dropped by 
performing this selection in the dereddened plane; they 
are fairly uniformly spread over the entire mapped region 
(and hence arc likely contaminants rather than YSOs). 
Just 3 of these are in the Gulf of Mexico region (sec R09) 



where we do not expect very much background contami- 
nation due to the very high reddening. Satisfied that this 
decision does not significantly affect our final sample of 
YSOs, we have chosen to keep our source selection in the 
observed plane. 

Once these background contaminants are removed 
from consideration, we are left with a population 
of 61,989 IRAC sources, dominated by objects with 
the apparent co l ors of stellar photospheres. As in 
IGutermuth eTaLl ([2008h . we have further selected YSO 
candidates in the [4.5] — [8] vs. [3.6] — [5.8] color-color 
diagram meeting the following criteria: 



[4.5] - [8] > 0.5 

& 

[3.6] - [5.8] > 0.35 

& 

I [3.6]-[5.8]<^x([4.5]-[8]-0.5) + 0.5 



(3) 



This 



This selection is indicated in Figure 31 panel b. 

leaves 1657 candidates. 

Continuing to follow IGutermuth et al.l (|2008f ). we used 
their selection criteria to identify objects from the 1657 
candidates whose IRAC fluxes may be contaminated by 
emission lines from shocks given their [3.6] — [4.5] and 
[4.5] — [5.8] colors by the equations: 



[3.6] - [4.5] > fii([4.5] - [5.8] - 0.3) + 0.8 



[4.5] 
[3.61 



& 

[5.8] < 0.85 

& 

[4.5] > 1.05 



(4) 



We found only 2 YSO candidates that matched these 
criteria (205608.3-1-433654.2 and 205702.1+433431.3). 
Their SEDs are compatible with their being real YSOs 
whose aperture photometry is affected by emission from 
a compact, circumstellar nebula, causing an excess at 
4.5 /xm relative to the two adjoining IRAC bands. More- 
over, these 2 objects arc located inside the "Gulf of Mex- 
ico" where the concentration of YSOs is the highest in the 
cloud, and where previous optical emission line surveys 
have found numerous HH regions. Given these consid- 
erations, we retain these objects in our YSO candidate 
list. 

Of the 1657 YSO candidates, 972 (59%) have a 2MASS 
counterpart. Out of the entire region, 131 (7.9%) have 
an optical counterpart; out of the region covered by the 
optical data, 10.4% have an optical counterpart. All 
but 63 of the YSO candidates are inside the region cov- 
ered by MIPS data, and 45% of the candidates have a 
MIPS 2Afim counterpart. We note that the percentage 
of 2MASS and optical YSO counterparts is less than for 
the whole IRAC 4-ba nd catalog (92% and 28% respec- 
tively, see section [53]) ; this is due to the fact that YSOs 
are mostly located inside highly reddened regions and 
are not detected in our shorter wavelengths. With this 
selection of ^^1600 YSOs, we have increased the num- 
ber of likely YSOs in the NANeb region by an order of 
magnitude. 

We compare the positions of both the YSO candidates 
and background contaminants in a 2MASS / IRAC color 
color diagram in panel d of Figure [4l A large fraction 
(90%) of sources flagged as galaxies or AGN fall in the 
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region occupied by reddened main sequence stars, and 
hence blueward in Kg — [4.5] color from the YSOs. Just 
5% of the YSO candidates are located in this same region, 
two-thirds of which are classified as intermediate between 
Class II and Class III (see SgS]). 

It is hard to estimate the contamination fraction of our 
YSO candidates without additional data (e.g. spectro- 
scopic or X-ray confirmation). However, we can make 
a few worst case estimates. Because our survey area is 
fairly large, it includes some sections which are outside 
the main star-forming region and have relatively few can- 
didate YSOs. The north-cast corner of our survey region 
(see Figure for example, offers one such relatively 
YSO-poor area. Specifically, we define a "field" region 
as being the area above a line from the mid-point of the 
East edge of Figure [3] to the mid-point of the North edge 
of the figure. If all the candidate YSOs in this region are 
contaminants, we derive a contaminant surface density 
of 32 objects per square degree. Assuming the AGN and 
other contaminants are uniformly distributed over our 
survey region, this yields an upper limit to the fraction 
of contaminants in our survey of 17%. 

3.2. NANeb age estimation 

We have used optical photometry in order to constrain 
the age of the NANeb complex. Figure [S] shows the loca- 
tion of YSOs with optical and 2MASS photometry in 
an optical color- magnitude diagram. Also shown are 
ISiess et all (j2000t) isochrones, where the color- Toff con- 
version has been tuned so that the 100 Myr isochrone 
follows the single-star seq uence in the Pleiades (jStauffeil 
[1991 [Jeffries et all [20071) . The number of YSOs in this 
diagram is limited by the fact that it requires optical 
counterparts. This means that most of sources from this 
diagram are located inside the Icss-extincted regions of 
the NANeb, on the edge of the Pelican Nebulae (80% 
are located inside reg ion of ^„ < 5, according to the 
iCambresv et al.l [20021 extinction map). The embedded 
Class I sources are not well-detected in the optical data, 
and our YSO selection criteria do not allow us to de- 
tect Class III sources, so this sample is mostly limited to 
Class II sources. 

Our age estimates are uncertain for all of the normally 
expected reasons (e.g. uncertainties in the isochrones 
and their transfer to the observational plane; the need 
for and imprecision of the reddening corrections; the ef- 
fect of spots and UV excesses on the estimated colors 
and effective temperatures of our target stars; binarity). 
However, it is apparent in Figure[5|that most of our YSO 
candidates (with optical counterparts) are younger than 
10 Myr according to the isochrone tracks. The median 
age is slightly older than 3 Myr, but the most embedded 
and probably youngest regions of NANeb are excluded 
from this optical CMD. The median age and the age dis- 
persion is compara ble to what found in NGC 2264 by 
IRebuU et al.l([200l . Since the reddening vector is essen- 
tially parallel to the isochrones, and YSOs with optical 
counterpart have low reddening, the median age will be 
minimally affected by reddening. 

3.3. Classes of YSO Candidates 

We classified this set of h ighly likely NANeb YSOs, 
according to the 'Lada' (|1987[ ) Class I/II/III system, up- 
dated bv lAndrd (,1994) where the "flat" class has been 



added between Class I and II. Classes are assigned based 
on the spectral index, defined by 

d logAF(A) 

a = 

d log A 

where F{X) is the flux at the wavelength A. For Class I, 
0.3< a; for flat sources, —0.3 < a < 0.3; for Class II, 
-1.6 < a < -0.3; and for Class III, a < -1.6. Note 
that because we are making a Spitzer-bascd selection of 
YSO candidates, by definition none of our candidates 
will have a photospheric slope, so our Class HI inventory 
is guaranteed to be incomplete. Further observations at 
other wavelengths {e.g. with Chandra) are required to 
find such objects. In the following, to avoid confusion, 
we call our YSO candidates with a < —1.6 Class II-III. 

We fit a over the 4 IRAC channel points from 3.6 to 
8 /im. From the total 1657 YSO candidates, we found 
that 1059 (64%) are Class II, 184 (11%) are flat and 198 
(12%) are Class I. We classified the remaining 216 (13%) 
YSOs as Class II-III because they have a between stellar 
photospheres {a = —2.7) and Class II (a = —1.6). A 
histogram of a for our YSO candidate sample is plotted 
in Figure [5] This figure also contains, for comparison, a 
histogram of the values of a for our entire NANeb catalog 
(mostly field stars), and that for the YSO candidates 
derived from the highly reliable Serpens catalog obtained 
bv lHarvev eFall ([2006, 2007b.a) (see discussion below). 

Plots similar to those appearing here (or in R09) can 
be found in most of the c2d series of papers on Ser- 
pens, Perseus, Ophiuchus, Chamaeleon II, and Lupus. 
We picked Serpens for our primary comparison here be- 
cause it is thought to be similar in age to the NANeb, it 
contains a deeply embedded cluster similar to the Gulf 
of Mexico, it is i n the Galact ic plane (^=32, 6=-|-5), 
and iHarvev et all (|2?M I2007bl [a[l performed a Spitzer- 
specific source selection. However, there are some impor- 
tant differences. The Spitzer maps of Serpens are only 
0.85 deg^, nearly 6 times smaller than our map. The 
selection m ethod for find i ng YSO can didates is much 
different in IHarvev et all l| 20061 l2007bl [a[). for example, 
requiring a MIPS-24 detection, and the differences are 
most strongly apparent in the selection of Class HI ob- 
jects. 

In order to fairly compare our data with the Serpens 
data, we used the highly reliable Serpens catalog (of all 
objects, not just their YSO candidates) provided by the 
c2d team as part of their final Legacy delivery (available 
on the SSC website) and applied our selection method 
to those data. From the sample of Serpens YSOs which 
pass our selection criteria, we calculated a in the same 
way that we did for the NANeb, from 3.6 to 8 fim. It is 
these values which appear in Figure [6| Our selection and 
classification scheme yield 52 Class I, 24 Flat disk, 189 
Class II and 65 Class II-III Serpens YSOs. The number 
of Class II-III o b jects is hard to compare to that found 
in IHarvev et al.[ ( 2006[ ) since they use the 24^m data in 
their classification scheme ; but the number of Class I, 
Flat and Class II are rough l y com parable: 30, 33, 163 
respectively in IHarvev et al.l (|2006( ). 

As a means to estimate a comparative evolutionary 
age for the NANeb stars, one can determine the ra- 
tio of the number of Class II ("older") to Class I or 
Flat sources ("younger"), and compare that ratio to 
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Fig. 4. — IRAC color-color and color magnitude diagrams used to reject background contaminants (panel a, b and c) and select YSOs 
candidates (panel h). An additional diagram shows the distribution of background contaminants and YSO candidates in a 2MASS-IRAC 
color-color diagram limited to targets with 2MASS counterparts (panel d). In all diagrams, blue triangles are sources flagged as PAH- 
emission sources (galaxies) located in the darkest gray area of panel a or b. Red squares are sources flagged as AGN, selected in the gray 
area of panel c. Finally, YSOs candidates are plotted with black cross symbols; they have been selected in the brightest gray area of panel 
b if not previously flagged as background contaminants. The reddening vectors correspond to an Ax of 10 in panels a, b and c, and Ak 
of 2 in panel d. We averaged the extinction law given in ll^'lahertv et all 1120071 ) for Serpens, Orion and Ophiuchus. 



that derived for other clusters, such as Serpens. For 
our entire NANeb YSO sample, we derive a ratio 
N(ClassII)/N(Classl4-Flat) 1059/382 = 2.78 ± 0.17. 
Using our similarly analysed data for Serpens, we derive 
this ratio as 2.49 ±0.33 (roughly comparable to what was 
derived by Harvey et al. 2006 - 2.6±0.4, where they used 
2MASS through 24 /im data to compute alpha). The ra- 
tios for NANeb and Serpens arc comparable, suggesting 
that the mean ages of the two YSO samples are indeed 

si milar. 

iHarvev et all ((20061 12007bllal l used a series of Spitzer- 
based selection cri teria to select a sampl e of YSOs; their 
final selection in IHarvev et al.l ()2007af ) primarily uses 
color-color and color-magnitude diagrams to assign a lik- 
lihood that a given object is extragalactic. The process 



starts with requiring detection in all 4 IRAC bands as 
well as MIPS-24. Only about half of our minimally con- 
taminated YSO sample has a MIPS-24 detection, so the 
Harvey et al. selection process by its very nature could 
only retrieve about half of our sample. However, '-^90% 
of our highly reliable sample with MIPS-24 detections is 
also retrieved by a Harvey et al. -style source selection. 

R09 contains a special discussion of the Gulf of Mex- 
ico, which is an interesting area containing hundreds of 
deeply embedded young stars, mostly associated with 
3 subclusters. We note here that our requirement for 
having all four bands of IRAC for our source selection 
omits many of the objects in this region, and this re- 
quirement will be relaxed in R09 to find cluster mem- 
bers. However, using our selection criteria and our 
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v-i 

Fig. 5. — y/y — / color magnitude diagram for our YSO candi- 
dates. The solid lines are, from right to left, tuned Siess models 
(see text) for 1 3 and lOMyr scaled to the NANeb distance. A 
reddening vector of Ay=5 is plotted. 




Fig. 6. — Normalized histograms of the distribution of spectral 
slope, a. The continuous line is the histogram for the YSO can- 
didates for the whole sample. All other IRAC sources from the 
NANeb are also indicated (in gray), as are the 3 30 YSO candi- 
dates we have found from the IHarvev et al.l 1)20061 ) highly reliable 
catalog (dashed line). 



best possible YSO sample (as defined above) results 
in a ratio of N(Class II)/N(Class I+Flat) that is sta- 
tistically significantly different inside the Gulf cluster 
than outside of it. The ratio, with its Poissonian er- 
ror, is 901/261=3.45 ±0.24 outside the Gulf of Mexico 
and 158/121=1. 31±0. 16 inside this region. Moreover, 
we find comparatively very few Class II-III objects in 
the Gulf of Mexico region. Both these findings indicate 
that the Gulf of Mexico cluster is, in evolutionary terms, 
the youngest region within the NAN complex. 

We explore now, how the NANe b reddening can chang e 
our YSO classification. Indeed iMuench et all l)2007l ). 
appendix A, demonstrated that for a deeply embedded 



cluster (At, ~40), the IRAC SED slope of a typical K6 
Class II member of IC 348 would be red dened into an ap- 
)arent flat spectrum. According to the iCambresv et al.l 
2003) extinction map, however, 97% of YSOs are lo- 
cate d in region whe r e the Ay is less than 10. Since 
the ICambresv et al.l ()2002[ ) map yields extinctions up 
to about 30, we assume that for extinctions < 10 the 
predictions should be reliable (i.e. not "saturated" be- 
cause the extinctions are so high that no background 
stars are present within the grid point). To test the ef- 
fect of reddening on our classification, we have dered- 
dened our SEDs prior to the spectral slope calculation, 
according to the Cambresy et al. map. The numbers 
of Class I, Flat disk and Class II given by the dered- 
dened SEDs - 186, 177, 1043 - are quite similar to those 
using the observed SEDs - 198, 184 and 1059, making 
the N(Class II)/N(Class I+Flat) ratio equal to 2.87±0.17 
instead of 2.78±0.17. This small differences does not 
change our main conclusions about the spatial distribu- 
tion of Class I and Class II. 

There is another means by which we can gauge 
the effect of ex t inctio n on our YSO classification. 
iGutermuth et al.l (|2008f l distinguished protostars and 
Class II by there [4.5] — [5.8] color excess. They chose 
this criterion specifically because it is relatively im- 
mune to extinction (because the interstellar extinction 
curve is relatively flat b etween these two wavelengths) . 
IGutermuth et al.l (|2008f l selected protostars as sources 
matching the following equation: 



[4.5] - [5.8] > 1 

OR (5) 
[4.5] ~ [5.8] > 0.7&[3.6] - [4.5] > 0.7 



The lGutermuth et al.l ()2008[ ) classification scheme does 
not include flat-spectrum sources as a separate class. Fig- 
ure [7] shows the location of Class I / Flat / Class II and 
Class II-III YSOs in the [4.5] - [8.0]/ [3.6] - [4.5] plan 
superimposed on the protostar area defined above. The 
location of our Class I and Class II sources agree well 
with the regions defined in this diagram: 95% of our 
Class I and 99% of our Class II sources would be classi- 
fied as protostars and Class II by the Gutermuth et al. 
criteria. Our fiat spectrum sources lie on the border be- 
tween Class I and C lass II, with roughly equ al numbers 
on either side of the IGutermuth et al.l ()2008f ) boundary. 



3.4. Spatial Distribution of YSO Candidates 

Now that we have a minimally contaminated sample 
of YSO candidates, we can look at their spatial distribu- 
tion and investigate the degree of clustering and the rel- 
ative spatial distribution of YSOs as a function of their 
class. We have computed a YSO density map of our 
obser ved region. We used a Kernel method (jSilvermanl 
|1986( ) which yields a smooth isodensity contour map from 
the projected position of objects. In each point of the 
space (a, S), the density is set, due to the contribution of 
all n points, by the kernel density estimator: 



1 

i=l 
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TABLE 3 
NANeb clusters 



[3.6] - [4.5] 

Fig. 7. — Spitzer color-color diagram. Open squares are Class II- 
III, plus symbols are Class II, green triangles are classified as Flat 
spectrum sources and circles are C lass I. The green area i s the 
protostar selection area defined by IGutermuth et al.l pOOSD . A 
reddening vector of Afc=10 is plotted. 

For the kernel, we adopted a 



(((5-*0' + (a-aO'cos2 6) 

2/^2 



where K is the kernel 
Gaussian shape : 

K{a,ai,S,6i) = -!- exp ■ 

where h is the smoothing parameter. Wc adopted h = 
0.05° (0.5 pc), which is approximately the size of the 
smallest group of YSOs in the NANeb discernible 'by 
eye'. 

Figure [8] shows contours of the density map we ob- 
tained, where both background contaminants and YSO 
candidates are plotted in two different panels. Using 
these density contours, we find first that the background 
contaminants (galaxies or AGN) are relatively evenly dis- 
tributed across the field of view. However, we notice a 
lack of background contaminant sources on the central 
dark cloud. The lack of contaminants in the high extinc- 
tion regions can be quantitatively explained simply due 
to the fact that our survey is magnitude-limited, and the 
extra extinction causes the observed fraction of galax- 
ies to drop out of our sample. We checked this assump- 
tion by creating an artificial map of randomly distributed 
background sou r ces re ddened by the extinction map of 
ICambresv et al.l (|2002[) . We then apphed a magnitude 
cut to this artificial sample, corresponding to the effec- 
tive faint limit imposed by our 4-band YSO selection 
criteria. The resulting spatial distribution of selected 
galaxies mimics well the observed galaxy distribution in 
Figure 7(top), confirming our assumption. 

Our YSO candidates are located primarily on the 
central dark cloud (Lynds 935), and are for the most 
part highly clustered. Further, we find that half of 
our YSO population is located in regions denser than 
1000 YS0s/dcg2 (~105 YSOs/pc^) and cover a relatively 
small fraction of the sky (^0.5 dcg^ compared to the ~9 
deg^ observed). 

From the distribution of YSOs, we distinguished by 
eye 8 main clusters (for a discussion of more quantitative 
means to identify clusters of YSOs from similar catalogs, 



id 


RA 


Dec 


Area 


YSOs 


< Ay > 


idc 




J(2000) 


J(2000) 


deg^ 


Count 


mag 




1 


20 56 17.4 


-1-43 38 18.8 


0.015 


48 


18.12 


1 


2 


20 57 07.1 


-1-43 48 21.8 


0.025 


120 


16.47 


2 


3 


20 58 19.0 


-1-43 53 32.5 


0.023 


74 


4.93 


3b 


4 


20 55 40.5 


-1-44 06 20.0 


0.025 


39 


5.68 




5 


20 50 54.9 


-1-44 24 18.1 


0.067 


121 


5.54 


5 


6 


20 53 35.2 


-1-44 27 57.4 


0.075 


130 


6.27 


6 


7 


20 59 14.2 


-1-44 16 41.3 


0.022 


22 


1.48 




8 


20 47 45.1 


-1-43 43 29.4 


0.017 


28 







The id, center coordinate and area of identified clusters arc indicated 
on the first 4 columns, followed by the number of YSO members and 
the average reddening from the extinction map of 
i2002l) . The last column is the cluster label from 
]2002D . if appropriate. 



Cambresv et a 



Cambresv et a~l 



see for example Jorgensen et al. 2008). We determined 
the center of each cluster as the local peak of the YSO 
density. We arbitrarily define the boundaries of these 
clusters as the YSO density contour at a level of 1/4 
of the maximum density peak. The location, the size, 
the number of stars and the mean Ay for each cluster 
is provided in Table [31 Five of these clusters have been 
previously identified bv ICambresv et ahl ()2002f ): they are 
the most prominent clusters located inside the Gulf of 
Mexico and in the central part of the NANeb. Figure [9] 
shows the location of both our clusters and clusters de- 
fined by Cam bresy et al. One c an see in this figure that 
cluster #4 of ICambresv et af] (|2G02D does not exist in 
our YSO candidate distribution. Clusters #7 and #8 of 
Cambresy appear in our YSO distribution, but are simi- 
lar to many other small clusters in the NANeb; they are 
not as densely populat ed as t he on es we have defined. 
Because the [Cam bres v et al.l (I2002D technique is sensi- 
tive to both WTTs and CTTs, whereas ours essentially 
excludes WTTs, it is reasonable that our two methods 
will yield somewh at different results. It i s likely that the 
clusters found bv ICambresv et all l|2002D but not by us 
may be the evolutionarily older clusters in the region. 

In the bottom panel of Figure [51 we plot the loca- 
tion of YSOs with different symbols corresponding to 
their class. To investigate in more detail a possible spa- 
tial segregation between Class I(-l-flat), II, and II-III, we 
used the YSO density contours to define regions within 
which we have computed the ratio of N(Class I -I- Flat) 
/ N(Class II) and N(Class II) / N(Class H-III). The be- 
havior of these ratios as a function of the density contour 
limit appears in Figure 1101 This Figure shows that in 
dense regions (i.e., more clustered regions), the propor- 
tion of Class II objects compared to the (presumably) 
younger Class I's is statistically lower than in the en- 
tire observed region by a factor of ~1.3. This is con- 
sistent with expectations - either due to high velocity 
stars moving away from their birth sites, or as a result 
of the overall expansion of young clusters as they age 
due to, for example, removal of their remaining molecu- 
lar cloud material by stellar winds. Moreover, the ratio 
N(Class II) /N(Class II-III) decreases from denser to less 
dense regions, as expected for the same reasons, by a 
factor of ~2. 

Figure [TOl also shows the fractional number of back- 
ground contaminants as compared to the fractional sur- 
face area inside each of the YSO density contours. The 
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Fig. 8. — TOP: the spatial distribution of background contaminants (red squares: AGN like sources, blue triangles: galaxies; see text 
and Figurc|4]l. BOTTOM: the spatial distribution of YSO candidates, where symbols denote the object Class (circles: ClassI or Flat, +: 
Class II, square: Class II-III). In both panels, wc plot contours of YSO density (see text for more details); the contour levels correspond 
to values plotted in Figure [10] from 200 to 2000 YSOs/deg-^. The contour which surrounds half of the total YSO population (~1000 
YSOs/deg^) is highlig hted in green and with a thicker contour line. The background gray area indicates the Ay >5 dark cloud from the 
ICambresv et al.l ^002 ) map. 
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Fig. 9. — Location of clusters shown on the distribution of the 
'^1600 YSO candidates (plus symbols). Solid red lines arc contours 
of clusters we adopted (see text for more details) . Our eight clusters 
are labelled with larg e, bold numbers and no brackets. The location 
of clusters found by ICambresv et all I I2002I) are also indicated by 
smaller bracketed labels. 



proportional number of background contaminants fol- 
lows very closely the proportional surface area. This 
suggests that sources flagged as contaminants are really 
background objects and not YSOs; otherwise one would 
expect an excess of sources flagged as background inside 
dense regions. 

Another way to characterize the degree of clustering, 
used by many authors, is the distribution of nearest 
neighbors. For each of our YSO candidates, we calcu- 
lated the distance to the 4th nearest YSO candidates; 
Figure [Tl] shows the histogram of this distance, displayed 



2000 1500 1000 500 

YSO density up 10 X YSOs/deg- 

FlG. 10. — Ratios of object classes as a function of YSO density. 
The X-axis corresponds to the YSO density contours plotted in 
Figure [8] from hight density to low density (corresponding to small 
and large surface surrounded by each contour on Figure 01 • Inside 
each contour, we calculated several fractional numbers: TOP: the 
fraction of YSO candidates (filled dot symbols linked by solid line) 
and background contaminants (filled triangle liked by dot-dashed 
line); the red dashed line shows the surface area surrounded by 
each contour as a fraction of the total covered surface (~9 deg^) 
MIDDLE: same, but for the fraction of (Class I-|-Flat) / Class II. 
BOTTOM: same, but for the fraction of Class II / Class II-III. 
Note that for the last value, X=0 corresponds physically to the 
boundary of our IRAC coverage. 

for Class I-I-Flat, Class II and Class II-III. Wc find a sta- 
tistically significant difference between the 3 histograms: 
using a Kolmogorov-Smirnov test, the probability that 
each of these distributions is drawn from the same par- 
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log io(d4th nearest) 

Fig. 11. — Histogram of the log of the distance of the 4th nearest 
YSO (in degrees) for each of our YSO candidates. The red dashed 
line is for Class I+Flat spectrum sources, the solid black line for 
Class II stars and the blue dash-dot histogram for Class II-III. The 
gray histogram represents the remaining population of non-cluster 
members. The Class I population, usually assumed to be younger, 
is significantly more clustered than Class II-III; this is consistent 
with expectations that, as stars dynamically evolve, they move 
away from their birthplace. 

ent distribution as either of the other two distributions is 
less than 0.02%. This is consistent with what we found 
above: the Class I population, usually assumed to be 
younger, is significantly more clustered than Class II-III. 

4. HH OBJECTS 

lOgura et al.l (|2002f l identified some HH objects using 
Ha over a small region of the complex near b right- 
rimmed cloud (BRC) 31. iBallv fc Reipurthl ([200l stud- 
ied the entire Pehcan Nebula using Ha, [N IT], and [S II], 
finding several new HH objects. 

The 3.5 to 8.5 iim wavelength range covered by 
the IRAC four channels is quite suitable to study 
deeply embedded young stellar outflows, since this 
range includes some of the brightest collisionally excited - 
pure rotational H2 emissio n lines (jWright et al.l 119961 : 
iNoriega-Crespo et al.ll2004al B) as well as other H2 vibra- 
tional lines (jSmith fc Rosen l2005l ) . At the longer wave- 
lengths, the MIPS channels sample a rich combination 
of atomic/ionic and molecular lines that include some of 
the best atomic species to study outflows, e.g. [Fe II] 
25.98 Mm, [O I] 63.18 Mm and [C II] 15 7.74 /xm (see e.g. 
iMoHnari et"alll2000l : iMorris et all 120041) that fall within 
the 24, 70 and 160 /xm bandpasses, respectively, plus of 
course continuum emission from cold dust. Wc discuss 
below our observations of one HH object discovered pre- 
viously, HH 555. 

HH 555 itself is not an embedded flow (like e.g. 
HH 211) since it is clearly detected at optical wavelengths 
(Ha and [S II] 6717/31 A) and its spectrum is consistent 
with that of a shock ex cited gas moving supersonically 
(jBallv fc Reipurthll2003l ). HH 555 belongs to that special 
class of irradiated lets, li ke HH 399 in the Trifid nebula 
(|Yusef-Zadeh et al] [2001 that are found in active star 
forming regions and are surrounded by a bath of UV 
ionizing photons from recently formed massive stars. 



The case for HH 555 is particularly interesting because 
the jet and counter-jet are bent westward as soon as 
they arise from their embedded source at the dense tip of 
the pillar, indicating the presence of stellar "sidewind" 
(|Masciadri fc Ragall200lD . The pillar itself is being pho- 
toevaporated l eading to the format ion of a two-shock 
structure (see iKaidic fc Ragal 120071 for details). The 
length of time that the jet will maintain its bipolar mor- 
phology is expected to depend on the strength of the 
impinging ionizing flux. The two-shock structure devel- 
oped by the interacting winds actually creates a shield 
that protects both the jet and the pillar, and there- 
fore, one would expect lower ionization on the West side 
of the outflow. Indeed, the Spitzer images combined 
with the ground based Ha emission (Figure [T2|) con- 
firms this stratification, if one assumes that the 8 yum 
and 24 /im emission arises mostly from H2 and [S II], 
as is the case in other jets. This also explains why 
the 24 /zm source lies in the middle of the 8 /im and 
24 M m jets, but is slightly offset with respect to Ha (see 
e.g. IKaidic fc Ra"gall2007l Fig. 3). 

Our 70 fim images do not show a clear signature of the 
jet, and this could be due to a combination of things: 
lack of sensitivity, artifacts along the scanning direction 
and physics ([O I] 63.18 /im is not present, neither is cold 
dust) . Careful analysis of the data do show the presence 
of the source at 70 at the tip of the pillar. The pil- 
lar itself, although faint, is distinguished and its emis- 
sion is likely to be due to [O I] photodissociated emis- 
sion. We measured the flux at 24 and 70 /im using a 
small aperture of 7" and 16" respectively, centered at 
20h51ml9.52s -h44d25m38.3s. We found 10±3 mJy at 
24 /im and 830±220 mJy at 70 /im. 

5. CONCLUSION 

We have combined mid-IR photometry from 
Spitzcr/IRAC with JHKs data from 2MASS to 
provide a sensitive photometric census of objects to- 
wards the North America and Pelican Nebulae. We 
have used IRAC color and magnitude diagnostics to 
identify more than 1600 sources with infrared excess 
characteristic of young stars surrounded by a disk. 

Our YSO candidates are located primarily on the cen- 
tral dark cloud (Lynds 935). We identified 8 main clus- 
ters, which together contain about a third of all the YSOs 
in the r egion; 5 of these cluste rs were previously recog- 
nized bv lCambresv et al](|2002l ). The sources wc identify 
as background galaxies or AGN are relatively uniformly 
spread across the observed region. These two properties 
suggest a low level of contamination in our YSO sample. 

We have assigned our YSOs classifications in the 
Class I/Flat /II/III system according to their infrared 
SED slope. The proportion of Class I+Flat compared 
to Class II is variable within the NANeb, with the dark 
clouds called "Gulf of Mexico" having a remarkably 
higher proportion of Class I's, suggesting these clusters 
have a younger age than the rest of NANeb. The "Gulf 
of Mexico" region will be discussed in more detail in R09 
where the MIPS data for NANeb will be presented. 

We compared optical photometry of a small sample of 
the YSO candidates to evolutionary models in order to 
provide an age estimate for the NANeb. We found that 
most of these YSOs app ear younger than 10 Myr accord- 
ing to the isochrones of iSiess et al.l (|2000f ). with the me- 
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Fig. 12. — Three color image of HH 555 arising from the tip of a dense pillar using Ha (blue), IRAC 8 ^m (green) and MIPS 24 ^m (red). 
The FOV is ~ 6' X 4', and North is up and East is left. The mid/far infrared bipolar jet (~ 1' in length) is centered on a 24 /im embedded 
source. The stratification in ionization along the jet from East-to-West (the gra dient in color) can be explained by the shield created during 
the two (stellar and photoevaporation) wind interaction IIKaidic &: Ragall2007l ). 



dian age being 3 Myr. However, the most embedded and 
probably youngest YSOs were not included in this sam- 
ple because our optical photometry is less sensitive than 
our IRAC observations and because the youngest stars 
are very red. However, despite all the caveats regarding 
this age estimation, it is the most robust published so 
far. 

Finally, thanks to our infrared images from 3.6 
to 70 fiui, we have provided new clues on the na- 
ture of HH 555, a pe culiar HH object reported in 
iBallv fc ReipurthI ()2OO30 . We detected a source at 70 /xm 
at the base of the bi-polar jets. Our data add support 
to the hypothesis that the bent shape of HH 555 is a re- 
sult of having the jet propagate in an environment with 
a stellar "sidewind" . 
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National Aeronautics and Space Administration and the 
National Science Foundation. These data were served by 
the NASA/IPAC Infrared Science Archive, which is op- 
erated by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the National 
Aeronautics and Space Administration. This research 
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the permission of these institutions. 

The research described in this paper was partially car- 
ried out at the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the National 
Aeronautics and Space Administration. 



REFERENCES 



Andre, P. 1994, in The Cold Universe, ed. T. Montmerle, C. J. 

Lada, I. F. Mirabel, & J. Tran Thanh van, 179 
Armond, A. C, Reipurth, B., & Vaz, L. P. 2003, Bulletin of the 

Astronomical Society of Brazil, 23, 92 
Bally, J. & Reipurth, B. 2003, AJ, 126, 893 
Bally, J. & Scoville, N. Z. 1980, ApJ, 239, 121 



Baraffe, I., Chabrier, G., AUard, F., & Hauschildt, P. H. 1998, 
A&A, 337, 403 

Cambresy, L., Beichman, C. A., Jarrett, T. H., & Cutri, R. M. 

2002, AJ, 123, 2559 
Comeron, F. & PasquaU, A. 2005, A&A, 430, 541 
Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 



15 



Flaherty, K. M., Pipher, J. L., Megeath, S. T., et al. 2007, ApJ, 
663, 1069 

Gutermuth, R. A., Myers, P. C, Megeath, S. T., et al. 2008, ApJ, 
674, 336 

Harvey, P., Merm, B., Huard, T. L., et al. 2007a, ApJ, 663, 1149 
Harvey, P. M., Chapman, N., Lai, S.-P., et al. 2006, ApJ, 644, 307 
Harvey, P. M., RebuU, L. M., Brooke, T., et al. 2007b, ApJ, 663, 
1139 

Herbig, G. H. 1958, ApJ, 128, 259 

Jeffries, R. D., Oliveira, J. M., Naylor, T., Mayne, N. J., & 

Littlefair, S. P. 2007, MNRAS, 376, 580 
Kajdic, P. & Raga, A. C. 2007, ApJ, 670, 1173 
Lada, C. J. 1987, in lAU Symp. 115: Star Forming Regions, cd. 

M. Pcimbert & J. Jugaku, 1-17 
Laugalys, V. & Straizys, V. 2002, Baltic Astronomy, 11, 205 
Luhman, K. L., Stauffer, J. R., Muench, A. A., ct al. 2003, ApJ, 

593, 1093 

Makovoz, D. & Marlcau, F. R. 2005, PASP, 117, 1113 
Masciadri, E. & Raga, A. C. 2001, AJ, 121, 408 
Molinari, S., Noriega-Crespo, A., Ceccarelli, C., et al. 2000, ApJ, 
538, 698 

Morris, P. W., Noriega-Crespo, A., Marlcau, F. R., ct al. 2004, 
ApJS, 154, 339 

Muench, A. A., Lada, C. J., Luhman, K. L., MuzeroUe, J., & 
Young, E. 2007, AJ, 134, 411 



Noriega-Crespo, A., Moro-Martin, A., Carey, S., ct al. 2004a, 
ApJS, 154, 402 

Noriega-Crespo, A., Morris, P., Marlcau, F. R., ct al. 2004b, 

ApJS, 154, 352 
Ogura, K., Sugitani, K., & Pickles, A. 2002, AJ, 123, 2597 
RebuU, L. M., Makidon, R. B., Strom, S. E., et al. 2002, AJ, 123, 

1528 

Siess, L., Dufour, E., & Forcstini, M. 2000, A&A, 358, 593 
Silverman, B. W. 1986, Density estimation for statistics and data 

analysis (Monographs on Statistics and Applied Probability, 

London: Chapman and Hall, 1986) 
Skrutskic, M. F., Cutri, R. M., Stiening, R., ct al. 2006, AJ, 131, 

1163 

Smith, M. D. & Rosen, A. 2005, MNRAS, 357, 1370 
Stauffer, J. R. 1996, in Astronomical Society of the Pacific 
Conference Series, Vol. 109, Cool Stars, Stellar Systems, and 
the Sun, ed. R. Pallavicini & A. K. Duprec, 305 
Werner, M. W., RoeUig, T. L., Low, F. J., ct al. 2004, ApJS, 154, 
1 

Wright, C. M., Drapatz, S., Timmermann, R., et al. 1996, A&A, 
315, L301 

Yusef-Zadeh, F., Biretta, J., & Wardlc, M. 2005, ApJ, 624, 246 



TABLE 4 
Table of YSO candidates^ 



lAU Designation 


RA (2000) 


Dec (2000) 


2MASS id 


B 


V 


/ 


J 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


SST205659. 3+434753.0 
SST205646.8+434609.0 


20 56 59.326 
20 56 46.765 


43 47 52.993 
43 46 08.980 


20565934+4347530 
20564675+4346089 


20.195 
19.981 


18.747 
18.023 


15.963 
15.385 


14.286 
13.226 


13.485 
12.180 


13.174 
11.783 


12.609 
11.653 


12.287 
11.485 


11.892 
11.244 


11.348 
10.392 


5.060 



^ Table will be presented in its entirety in the electronic version. Headings are shown here for guidance regarding its content. 



